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Abstract 

We comment on the method of Dreisigmeyer and Young [D. W. Dreisig- 
meyer and P. M. Young, J. Phys. A 36, 8297, (2003)] to model nonconserva- 
tive systems with fractional derivatives. It was previously hoped that using 
fractional derivatives in an action would allow us to derive a single retarded 
equation of motion using a variational principle. It is proven that, under cer- 
tain reasonable assumptions, the method of Dreisigmeyer and Young fails. 
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1 Introduction 



In 1931 Bauer proved the following corollary 0: 'The equations of motion of a 
dissipative linear dynamical system with constant coefficients are not given by a 
variational principle' . There are a few methods that may allow us to get around 
Bauer's corollary. For example, we could allow additional equations of motion to 
result. Bateman used this technique in 0. If we use the Lagrangian 

C 

L = mxy + — (xy — xy) — mu 2 xy , (1) 

where C is a constant, we would have the following equations of motion 

mi + Cx + mui 2 x = (2) 
my — Cy + mu 2 y = . (3) 

Bateman's method uses the loophole that Bauer's proof assumed that no additional 
equations arise. 

Riewe pointed out that Bauer's proof also implicitly assumes that all of the 
derivatives are integer ordered @. This has led to attempts to use fractional 
derivatives in the actions to model nonconservative systems 0|5j|5|. Here we will 
close this second loophole by extending Bauer's corollary to include fractional 
derivatives. 

Our paper is organized as follows. In section E] we review the background 
material needed for our result. The extension of Bauer's corollary is proved in 
sectional A brief discussion follows in section^ 



2 Background material 

Here we develop the relevant mathematics for our proof. A fuller discussion of 
this material can be found in [4]. Fractional derivatives can be defined using the 
theory of distributions. First, define the generalized functions 

$+( t ) = J T(a) 1 1 > U (4) 
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and 

*«-(*) = { f HI (5) 

where T(a) is the gamma function. The left fractional derivatives (LFD) of a 
function q(t) is given by 

a Bf[q] := * (6) 

r(-«) A 

where we set q(t) = for t < a. When a = n, n an integer, © becomes 

oPt Wi = V n q (7) 

where V is the generalized derivative. Right fractional derivatives (RFDs) are 
defined similarly 

Pi[q] := $Z a (t) * q(t) (8) 

r(-«) y t 

where now q(t) = for t > b. Instead of ©, we have 

t T> n b [q] = (-l) n V n q. (9) 

In H the actions were treated as Volterra series. The Volterra series are a 
generalization to functionals of the power series of a function. For a functional 
V[q), define the symmetric kernels 

mn,...,r n ) := - f V[q] . (10) 

Now introduce the notation 

K^*q n := ! ■■■! Ki s \r l ,...,r n )q(T n )---q(r 1 )dr n ---dr 1 . (11) 

Then V[g] can be written as 

Vfe] = ^ 1 ^ s) ^ n (12) 

n=l 
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(si 

where we set Kq = 0. 

The $^ (£) are now treated as kernels in a Volterra series. We can then take the 
functional derivative of the series to derive our equations of motion. An example 
should make this clearer. We will restrict our attention to the action 

V[q] = \K 2 *q 2 (13) 

where K 2 (t, r) in (fT3l is an arbitrary kernel, i.e., not necessarily symmetric as in 
(fTOb . (Equation (fT3l would be sufficiently general to handle the nonconservative 
harmonic oscillator.) Now let K 2 (t, r) be given by 

K 2 (t,r) := m$+ 2 (t-r) + mC$l 7 (t - r) + muo 2 $+(t - r) (14) 

where < 7 < 2 and C is a constant. So (fT3l becomes 

V[q) = y JJ^t 2 (t-r) + C^(t-r)+uj 2 ^(t-r)} x (15) 
q{r)q{t)drdt . 
The functional derivative of (1731 is 
*V[g] 



Mp) 



777 / 

- y - r) + C$i 7 (p - r) + w a $+(p - r)] g(r)dr (16) 



retarded 



HI f 

+ y jf [*la(* - P) + CSt^t - P) + - P)] 9 



(t)dt 



advanced 

If we require the advanced and retarded parts of ([TBI) to vanish separately, we have 

m [$± 2 {t) + Cfct^r) + u; 2 $+(r)] * qr(r) = (retarded) (17) 
m [$I a (t) + C$I 7 (t) + u; 2 $ (t)] * g(t) = (advanced) (18) 

where in d 1 8b we used the fact that — r) = $~(r — t). This is the method 
presented in [4| for deriving the equations of motion. 



3 The result 



In section El we reviewed the procedure Dreisigmeyer and Young proposed in 
for deriving a system's equations of motion. From ( fT7t and (fTSl) we see that two 
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equations are actually derived: an advanced one and a retarded one. So this is, 
effectively, a generalization of Bateman's method (see - ©). That is, extra 
equations of motion are allowed to result from the action's variation. 

We desire to have a single, retarded equation of motion to result from a 
variational principle. From (TBI) we see that the derivative operators are always 
contained in the /^(ti,^) kernel. Perhaps it is possible to use some kernel 
other than the $J(ti — r 2 ) to have a fractional derivative arise from an action's 
variation? The following theorem shows that this is not possible within the 
formalism presented in @. 

Theorem 3.1 There does not exist a K(t\ , T2), t\,t% G K, such that the variation 
of the quantity 

V[q] = J K^rMnhiT^dndTz (19) 

will result in a D" [q] for a ^ 2n, n an integer. 
PROOF . The variation ofV[q] is given by 

|S = [K{p,t)+K{t,p)]*q{t). (20) 

We will assume that 

[K{p,t) + K(t,p)]*q(t) = *± a (j>-t)*q(t) (21) 

and arrive at a contradiction. We require that A211 holds for every q(t). Then we 
must have 

[K(p,t)+K(t,p)) = ®± a (p-t). (22) 

Interchanging p and t in X22[ gives us 

[K(p,t)+K(t,p)] = (23) 

Hence, unless $^ a (p — t) is symmetric in p and t, H22\) and (123\) cannot both hold. 
That is, unless a = In, n an integer, there does not exist a K(t\, t 2 ), T\, r 2 e M, 
such that XUl holds. ■ 
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Theorem 13 . 1 1 show s that, in general, the fractional mechanics formalisms pre- 
sented in 00 cannot derive a single, retarded equation of motion. In order to 
overcome this difficulty, Riewe suggested approximating RFDs with LFDs JUE1. 
Dreisigmeyer and Young showed in |4| that this is not a sound idea and, instead, 
allowed for an extra, advanced equation of motion. This latter technique is not, 
itself, entirely satisfactory. 



4 Discussion 

Theorem l3.1l shows that some revision of our concept of an action may be in order 
if we desire a variational principle to work for nonconservative systems. How 
could we derive a single, retarded equation of motion for systems? Our result 
holds even if K 2 (ti, r 2 ) is allowed to be complex. We would also require that 
(7(7-1) = q{r 2 ) for r\ = r 2 in (THJ1). That is, we do not want to employ Bateman's 
method, as was done in 0. 

One possible method proposed by Tonti [7J (see also JT)) is to use the convo- 
lution product in our Lagrangians. This leads to actions of the form 

V[q] = J K(t-n- T 2 )q(T 1 )q(T2)dT 2 dT 1 dt. (24) 

This method does allow the derivation of a single retarded equation of motion for, 
e.g., the driven nonconservative harmonic oscillator. Unfortunately, it does not 
seem possible to naturally generalize Tonti's method to higher ordered potentials. 
That is, using quantities like 

/ Kit - Ti - ... - T n )q(n) ■ ■ ■ q(r n )dr n ■ ■ ■ dndt (25) 

in the action will not lead to the correct form for the potential energy terms when 
n > 2. This situation should be contrasted with that in \4\. There terms like 

/ $0 ( T i ~ T 2) • • • $0 ( r "-i - r n)q(n) ■ ■ ■ q(r n )dr n ■ ■ ■ qn (26) 

were able to treat the potential energy terms correctly. However, as Theorem 
demonstrates, the formalism in is unable to deal correctly with the fractional 
derivative terms. 

An interesting feature of d2"4l versus (THJ1) is the presence of t in the action along 
with Ti and r 2 . Our stated goal is to derive purely retarded equations of motion 
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using a variational principle. However, to achieve this we need to find the correct 
kernels for our Volterra series action. Theorem l3.1l tells us that kernels of the form 
^2(71, T2) are not sufficient for our purposes. Equation d2"4l suggests that we look 
instead at the kernels K(t, n, . . . , T n ) for our actions. 

It is desirable to be able to use the same type of kernel for the fractional 
derivatives as well as the higher ordered potential terms. This assumption is 
necessary so that the same perturbation of q(r) can be used in the kinetic and 
potential energy terms of the action. It allows us to reject using kernels of the 
form K(t — n — ■ ■ ■ — r n ), as in d25l . Also, it prevents us from using terms 
like d2"4l for the kinetic energy and terms like d2~6l ) for the potential energy, within 
the Lagrangian formalism. Finding the correct kernels K(t, n, . . . ,r n ) for our 
Volterra series actions is a line of research that we are actively pursuing at this 
time. 
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